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Abstract 

We have compared the ground-state energy of several observed or proposed 
"2^/2 x 2^2 oxygen (O) ordered superstructures " (from now on HS), with 
those of "chain superstructures" (CS) (in which the O atoms of the basal plane 
are ordered in chains), for different compositions x in YE^CusOg+z- The 
model Hamiltonian contains i) the Madelung energy, ii) a term linear in the 
difference between Cu and O hole occupancies which controls charge transfer, 
and hi) covalency effects based on known results for t— J models in one and two 
dimensions. The optimum distribution of charge is determined minimizing 
the total energy, and depends on two parameters which are determined from 
known results for x = 1 and x = 0.5. We obtain that on the O lean side, 
only CS are stable, while for x = 7/8, a HS with regularly spaced O vacancies 
added to the x = 1 structure is more stable than the corresponding CS for 
the same x. We find that the detailed positions of the atoms in the structure, 
and long-range Coulomb interactions, are crucial for the electronic structure, 
the mechanism of charge transfer, the stability of the different phases, and 
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the possibility of phase separation. 
PACS numbers: 74.72.Bk, 72.15.-v, 61.66.-f, 74.10. +v 



I. INTRODUCTION. 

The ordering of oxygen (O) atoms in the basal plane of YP^CusOe+x, and its relation 
with electronic properties, particularly the superconducting critical temperature T c , has 
been a subject of great interest. Overviews are contained in Ref. Jl|. By the end of the last 
decade, for x near 0.5 evidence or ordering in chains was found, while for x near or 1, 
HS were reported 0-f|||, the charge distribution in superconducting planes || and in the 
Cu ions of the basal planes PJlOl was determined, and an explanation of these experiments 



was presented using a lattice-gas model based on Coulomb repulsions, and the appropriate 
extension to this system of the three-band Hubbard model JIT]]. Optical P,l2| and 



nuclear quadrupole resonance (NQR) p3| , |T3|| experiments give strong evidence that two-fold 
coordinated Cu atoms are in an oxidation state Cu + , while three- and four-fold coordinated 
Cu atoms are Cu +2 . Qualitatively, this fact is easy to understand in terms of |T^]: the 
energy necessary to add a d hole to an n-fold coordinated Cu + (surrounded by n 0~ 2 ions) is 
ecu — ZnUpd, where U p d is the Cu-0 hole-hole repulsion. Provided that the hole Fermi energy 
of the superconducting Cu0 2 planes lies between e C u — 4t/ p d and e Cu — 6U p d, only two-fold 
coordinated Cu atoms remain Cu + , and the remaining holes go to the Cu02 planes (and to 
O atoms of the chains). This mechanism not only provides the holes for superconductivity 
of the CS, but also reduces their energy |14[]. For x ~ 0.5, annealing at room temperature 
produces an increase of T c due to an increase in the amount of two-fold coordinated Cu ions 

EMU- 

In spite of this qualitative understanding of the relation between O ordering and charge 
transfer, and the success of lattice-gas models in explaining the basic features of the ther- 



modynamics of O ordering at temperatures above room temperatures pT| , |T9| -f2~2]j , a fully 



self-consistent theory of the atomic and electronic structure of YBa 2 Cu30 6+:! ; is still lacking, 
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and several controversial issues remain. Schleger et al. showed that it is necessary to add 
electronic degrees of freedom to lattice-gas models in order to explain the observed dx/dfi, 
where /i is the O chemical potential ||23|| . The strong-coupling approaches to the electronic 



structure are able to treat adequately the on-site Cu Coulomb repulsion Ud, and to explain 
the observed dependences of the hole count in the planes tin (from which T c can be inferred 
), and the amount of Cu + [|li p2| , p4| -[26|] . However, these agreements are obtained within 



a region of parameters of the model, which is not fully justified, and Coulomb repulsions 
beyond ~ 2A are neglected. On the other hand, ab initio calculations, as a consequence 
of the large value of Ud and the neglect of correlations, fail to describe the semiconducting 
phases (they predict a metal) and the neglected correlation energies are much larger than 
the lattice-gas model parameters which determine the structure IJ. 

A controversial issue is the stability of HS. While it is clear that at room temperatures, the 
CS disappear in the semiconducting region the experimental and theoretical situation 
does not allow at present to disclose unambiguously the nature of the ground state. The 
HS observed by transmission electron microscopy |@||, were proposed to be metastable 
0. Synchrotron radiation studies for x ~ 0.2 provided strong evidence of the presence of 
~ 0.2% of a parasitic phase BaCu304 [fZTJ], which is able to explain the x-ray diffraction 
pattern observed for x ~ 3/8, ascribed previously to O ordering ||28|| . This pattern is also 



fully compatible with the superstructure of minimum Coulomb repulsion between equally 
charged O atoms, among all those with unit cell 2^/2 x 2v^2 |28| (see Fig. 1). On the other 
hand, neutron-scattering experiments for x ~ 3/8 pSfl , which are more suitable to study O 
ordering due to the comparatively larger O cross section ||30|| , are so far only explained in 
terms of an O ordered superstructure of unit cell 2\/2 x 4 a/2 []30| , |31]1 . On the theoretical side, 



there is no general physical argument from which one can discard or confirm HS p|, |30| , |32| , |33 



except under too restrictive hypothesis |55 |, as we will show in Section III. The model which 
obtained these superstructures |TT| , |3l| is based on Coulomb repulsions between any two O 
ions, with metallic and dielectric screening and a parameter AE, which favors CS, to take 
into account the relation between charge transfer and Cu coordination mentioned above 



|I|j3^|. AE has been calculated using an extended Hubbard model, but the result is too 
sensitive to the parameters of this model, not known with enough precision. For large but 
reasonable AE, the ground state of the structural model is a CS for any O content x p|.|32|. 

Because of the very low or vanishing density of carriers, the Madelung energy is a funda- 
mental ingredient in the physics of YE^CusOg+x P5|-pS[ and other [38| high-T c systems. In 
YBa2Cu306+x, 0-0 repulsions at distances of at least 8A are necessary to explain observed 
split diffuse diffraction peaks f39 |, and 0-0 repulsions at distances ~ 27 A are required to 
stabilize several of the observed CS f32[. Also, the model of Ref. fllTf , based on Coulomb 
repulsions, seems to be able to explain qualitatively the experimental data gathered on 
YBa 2 Cu30 6+:! ; [|IJ, including the structure transformation kinetics if supplemented by long- 
range elastic energies |f40f . Here we generalize previous works considering the Madelung 
energy, allowing for the possibility of charge transfer and finding the optimum distribution 
of charges by minimization of the total energy. Effects of covalency are also included. Ne- 
glecting the latter, our approach is equivalent to the exact solution, in the limit of zero 
hopping and infinite on-site Coulomb repulsions, of the appropriate model of the Hubbard 



is 



type for the system. The effect of correlations, essential for Cu-0 charge transfer [TJ 
adequately retained. The main shortcomings of the approach are the neglect of core-core 
repulsive energy (which amounts to 10% of the Madelung energy and stabilizes the lattice) 
and screening effects. 

The theoretical treatment is presented in Section II. Section III contains the results and 
Section IV is a discussion. 



II. THE MODEL 

The ground-state energy of the system per Y atom is described as: 

E = E Mad + E A + 2E p i + E ch . (1) 

Euad is the Madelung energy as a function of Cu and O charges. is the energy required by 
the charge-transfer process Cu + + O — > Cu +2 + 0~ 2 in absence of interatomic repulsions. 
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E p i describes the kinetic and magnetic energy gain due to covalency in the superconducting 
Cu0 2 planes (assumed equal for both planes of the unit cell), and E ch is the corresponding 
term for the CuC^+z subsystem containing the basal plane. 



The Madelung energy can be written as pEIfl 



Euad = Pi^i'i Pi = a ijZj, (2) 

* j 

where e is the elementary charge, N the number of Y atoms in the supercell, the charge 
of the i th atom in the supercell, and $ the electrostatic potential at the position of this 
atom. Because of charge neutrality J2i%i = 0, the geometrical coefficients aty can be shifted 
by an arbitrary constant. We have chosen it in such a way that an = for all i. Then for 

1 ^ , 1 1 



a 



E( )• ( 3 ) 

r j,0 — r i,0 T ^o r i- T ~ r i,0 r i,T — r i,0 

Here T labels the translation vectors which map the superlattice onto itself and r^. T = 
T + rk,o is the position of the k th atom of that supercell obtained from the one lying at the 
origin by a translation T. The sum over T is evaluated by the Ewald's method |[42|| . For 
simplicity we have assumed that the lattice parameters a = b (taking the average between 
them), and we have taken the positions of the atoms from three different O contents:x = 7 
f43|| , x = 0.45 |Q and x = [^4| |. We also assume that all Y ions are Y +3 , all Ba ions are 



Ba +2 , all Cu ions of the superconducting planes (Cu(2) in the notation of Refs. |43|j44]| ) have 
the same charge, and all O atoms of these planes (denoted 0(2), 0(3) f4~3| , |4"4|| ) related by 
translations of the primitive lxl cell have the same charge. This simplifies considerably the 
problem, allowing us to express part of the sums in i and j of Eqs. (2) and (3), in terms of 
the coefficients a°- of the primitive unit cell, reducing appreciably the number of ct's which 
should be calculated. Neglecting an unimportant constant we can write: 

E A = 3A(Z Cu -2), (4) 

where Zq u is the average charge of all Cu atoms. In principle, A is the difference between 
the ionization potential of Cu + and the (negative) electron affinity of 0~. However, it should 



also contain information of steric effects (short-range repulsions), the energy gain of a local 
Zhang- Rice singlet |3(J (not included in E pt + E ch ) and the kinetic energy of the Cu hole. 
We keep A as a parameter, assumed independent of x and determined in such a way that 
the charge distribution in YBa2Cu306.5 agrees with experiment. 

The number of added holes in one of the two superconducting planes per unit cell is 
h = 2 + Zcu2 + Z02 + Z03, where Zc U 2, Z02, Z03 are the average charges of the Cu(2), 0(2) 
and 0(3) atoms of that plane. For the kinetic energy as a function of h, we take the form 
established from a high-temperature expansion of the t — J model, with t=OA eV, J=0.1 
eV [45|, slightly generalized to give the correct magnetic energy for h = 0: 

E pl = -[eh + (1 - h){h + 0.09192)eV]. (5) 

e represents the difference between the energy gain in forming a localized Zhang- Rice singlet 
in the planes with respect to forming it in perfect CUO3 linear systems (involving the Cu(l), 
0(4) atoms of the basal plane, and two 0(1) apical atoms p3| , ^4]] ). It can also contain infor- 
mation of different steric effects in chains and planes. We keep e as a parameter determined 
from the experimental charge distribution for x = 1. Actually, a realistic one-band model 
for the cuprates contains also hoppings beyond nearest neighbors and three-site terms which 
determine the shape of the Fermi surface and are critical for the superconductivity H46,47j, 
but we expect that these terms do not affect E p i very much. 

is obtained by fitting exact results for the one-dimensional t — J model. The values 
t = 0.85 eV, J = 0.2 eV were determined from a low-energy reduction procedure which leads 
to excellent results for the optical conductivity of the chains [i"8[|. The result is [ji8,33[|: 

E ch = [(J - 2t) sin(7r(l - h c )) - J(0.69 + 0.41/i c )(l - H c )} Vl - Jy 2 , (6) 

where y\ is the concentration of 0(4) atoms belonging to perfect chains, is the concen- 
tration of isolated 0(4)~ 2 ions between two Cu(l) +2 ions, and h c is the number of holes per 
Cu, added to (Cu0 3 )~ 4 perfect chains. This expression does not take into account doped 
extremely short chains or the correction corresponding to chains of intermediate length. 



The total energy E is minimized with respect to the charges in the subsystem Cu02+ X - 
(containing Cu(l), 0(1) and 0(4) atoms), and the average Cu(2), 0(2), 0(3) charges on 
both superconducting planes. Cu charges are allowed to vary between 1 and 2, and O charges 
between -2 and -1. This corresponds to the limit of very large on-site Coulomb repulsions. 



III. RESULTS. 

First, we have applied our approach to stoichiometric YBa2Cu307, a metallic state with 
a reasonably well-known, non trivial charge distribution. This is difficult to obtain from the 
Madelung energy, because this concept was developed for insulators. In fact, if we neglect 
E p iandE c h, we obtain that all Cu ions are Cu +2 , and all O ions are 0~ 2 except the chain 
0(4) atoms, the oxidation state of which is -1. In other words all holes go to the one- 
dimensional Cu(l)-0(4) chains which are not expected to conduct (due to defects or Peierls 
distortions) and the system is insulating. This is not bad as a first approximation. It is the 
best description of the observed charge distribution in terms of integer charges. Including 
all terms in the energy, with A < 46 eV and a reasonable e = 2 eV, we obtain that 60% 
of the holes enter CuO chains and 20% are in each of the superconducting Cu0 2 planes, 



in agreement with optical conductivity measurements |49(], and estimations based on bond 
valence sums and other experiments |50| . The resulting charges, Madelung potentials $ and 
different contributions to the energy are included in Table I. 

After checking that in general this gives lower energy, and to simplify the algorithm, we 
have constrained the minimization procedure distributing the holes on both superconducting 
planes in equal amounts between the 0(2) atoms of one plane and the 0(3) of the other, and 
have kept -2 the charge of the apical 0(1) ions. The results depend significantly, but not 
dramatically on e: for e=0, the amount of holes in each plane is reduced to 14%. Instead, 
the positions of the Ba and apical atoms are crucial. If the atomic positions are taken as 
those for x = 0, for which the 0(1) atoms are nearer and the Ba atoms more distant from 
the basal plane, we obtain a hole doping h of only 0.03 in each Cu0 2 plane, while for the 
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positions corresponding to x — 0.45, the resulting doping is 0.14. 

To establish bounds on A, we have calculated next the energy and charge distribution of 
both superstructures shown in Fig. 1 for x — 1/2. The HS is the one which minimizes the 
Coulomb energy when all atoms related by symmetry operations of the tetragonal primitive 
unit cell have the same charge [|I^,^T|,^4]] . However, the ground state is the CS, the two-fold 
(four-fold) coordinated Cu ions are mainly Cu + (Cu +2 ) [pj |i2t[l0i|l3| , and about 0.1 holes 



per Cu go to the superconducting Cu02 planes |||50||, in agreement with theory [0,^2]. To 
satisfy this charge distribution, our model has to satisfy several constraints. One of them is 
A > A' min = e(Pop — Pcu2) — dEpi/dh, where (3o P is the lowest $ of the occupied O atoms 
of the planes, and f3cu2 is the potential at the two-fold coordinated chain Cu(l) atoms. If 
this constraint is not satisfied, the holes of the planes go to the two-fold coordinated Cu 
atoms and the system would be insulating. In Table II we give the resulting charges, /3j and 
energies of the CS for A = 31 eV, slightly above A' min . For larger values of A (not too large 
to avoid that four- fold coordinated Cu +2 becomes Cu + ), the only change is that E and 
decrease proportionally to A/2. The resulting amount of holes in each Cu02 plane (0.093 
per Cu) is in very good agreement with experiment. The /3j at inequivalent 0(1) atoms 
are surprisingly similar and a little bit smaller than the $ of the O atoms of the Cu02 
planes. Holes prefer the latter because of the positive value of e. However a more realistic 
description should allow that a small amount of holes enter apical 0(1) atoms, particularly 
for small e. 

The experimental evidence indicates that three-fold coordinated Cu(l) atoms are Cu +2 , 
in particular optical experiments on quenched samples 0,0, and NQR experiments in which 



Y is replaced by larger ions fl3| . This is also suggested by the theoretical studies ||T|,|22|j3l 
and the argument on the interplay between charge transfer and structure [0] presented in 
Section I. An upper bound on A can be obtained requiring that in the HS for x = 1/2, all 
Cu ions remain Cu +2 and no holes are transfered to the Cu02 planes, so that the system 
remains semiconducting. This implies A < A Max = e(Pop ~ Pcu(i)) ~ dE p i/dh\ h=0 . We 
obtain A^ai = 38.29 eV. It is interesting to note that in order for the present work to be 
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consistent with the estimations of Ohta et al. for the charge transfer gap in several high- 
T c superconductors |52| one should take A ~ 3.5 x 10.9 eV = 38.15 eV. The results for 
any A < Amo,x are presented in Table III. Comparison with the energy of the CS (Table 
II), establishes a better lower bound for A: in order for the CS to be the ground state, 
A > A min = 31.34 eV. 

In the following, we assume that e = 2 eV and A min < A < Auax for all x. Using these 
criteria we derive conclusions regarding the stability of HS in comparison with CS. Let us 



begin with x = 1/8, calculated with the atomic positions for x = [ 44| . For CS, all four- fold 
coordinated Cu(l) ions remain Cu +2 , the two- fold coordinated ones are Cu + , and most of 
the holes apported by the neutral O atoms entering the x = structure to form that of 
x = 1/8, remain in their neighborhood: one hole is transfered to a nearest neighbor Cu + 
and only 0.09 additional holes per supercell are distributed in the planes. In other words the 
charge of the chain 0(4) atoms is -1.090 and the doping of each superconducting plane is 
h = 0.0056. The potential at the four-fold coordinated Cu(l) atom in the supercell is -24.55 
V, while those at two-fold coordinated Cu(l) atoms vary between -13.24 V to -12.36 V, with 
increasing distance to the Cu(l)-0(4) chain. The difference of more than 11 V is not taken 
into account in Hubbard-type models which do not include a large nearest-neighbor Cu-0 
repulsion U p d in an appropriate way [||.|I]].[33|]. A at the chain 0(4) atoms is 13.12 V. The 



j3i at the Cu02 planes and apical O atoms have similar values as those reported in Tables I 
and II. 

In the HS for x = 1/8, all O ions are O -2 and all three-fold coordinated Cu(l) ions are 
Cu +2 . fa at these atoms is -23.82 V, while at the remaining, two-fold coordinated Cu(l) + 
ions, it is ~ —12V. At the 0(4) ions = —18.30 V and other are similar as those in Table 
III. For A = A min = 31.34 eV, the total energy of the CS is -296.33 eV, slightly less than 
that of the HS (-296.07 eV). Although the Madelung energy of the latter is less (-272.57 eV 
in comparison with -268.63 eV of the CS), the CS has lower energy because it ionizes half 
of the Cu(l) + ions of the x = structure, and thus, pays less A. For larger values of A, the 
difference increases linearly with A/8. Thus, the CS is the ground state for x = 1/8. 
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Under some general conditions, using a multiband Hubbard model including U p d, one of 
us has shown that HS-type superstructures have less energy in the semiconducting phase 
if no holes enter apical 0(1) atoms |[33|| . What is the reason of the discrepancy with the 
present result? On one hand, the effect of repulsions beyond nearest neighbors is important. 
For example, in usual Hubbard-type models, the energy necessary to add a hole in a chain 
0(4) atom with both nearest neighbors being Cu +2 is e p + 2U p d independently of the rest of 
the electronic and atomic structure. However, as explained above, this energy is 13.12 eV 
for the CS, and 5.28 eV larger for the HS. On the other hand, to obtain the present charge 
distribution of the CS with the model of Ref. f33fl , e p < €d + 2U p d is required, contrary to what 
is expected in Cu0 2 planes p3| , |54]| , and one of the hypothesis of Ref. As stated clearly in 



Ref. |33[], that calculation was aimed to discuss the effects of covalency neglecting repulsions 



beyond nearest-neighbor Cu-0 ones and assuming that these were small. However, we find 
that longer range repulsions are essential. 

Next, we analyze the superstructures corresponding to x = 3/8 with atomic positions 
taken from data for x = 0.45. The general trends of the charge distribution and potentials 
are similar to those of x = 1/8. Two-fold coordinated Cu(l) ions remain Cu + , while higher 
coordinated Cu ions have an oxidation state Cu +2 . For CS, the charge of chain 0(4) atoms 
is -1.37 and the doping per Cu of the superconducting planes is h = 0.07. The potential at 
the sites of two 0(4) atoms of the unit cell is 15.01 V and 16.12 V at the remaining 0(4) 
sites. The $ at four-fold coordinated Cu(l) sites are near -25 eV, and those at two-fold 
coordinated Cu(l) sites vary between -11.45 V and -13.17 V. For the HS, all O ions are O -2 , 
the (3i at 0(1) atoms is near 20.8 V, and those at four-fold (two-fold) coordinated Cu(l) 
atoms are near -22 V (-10.35 V). We should note that keeping this chargedistribution, there 
are at least two superstructures with less energy [BTfj , one of them which provides the best 



fit of the neutron scattering data The difference between the Madelung energy of 

the superstructure of Fig. 1 and the lowest lying of the above mentioned superstructures 
is 0.17 eV Including this correction, the energy of the lowest lying superstructure of 
unit cell multiple of 2y/2 x 2^2 for A = A min = 31.34 eV becomes -294.62 eV, only slightly 
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smaller than the energy for the CS (-294.56 eV). Since increasing A favors the CS by a term 
proportional to 3A/8, there is a crossing already at A c = 31.50 eV and for A c < A < A Max 
= 38.29 eV, the CS has lower energy. If the x = 3/8 superstructures are analyzed with 
atomic positions corresponding to x = instead of x = 0.45, the same trends are observed. 
A c increases slightly to 31.85 eV. The doping of the planes for CS is reduced to h = 0.02. The 
energy decreases in ~ 0.6 eV for both structures. This difference becomes significant when 
one considers the possibility of phase separation: the structures with x = 3/8 calculated 
with the positions for x = (x = 0.45) [44| are stable (unstable) against phase separation 



into phases with x = 1/8 (calculated with the positions for x = 0) and x = 1/2 (calculated 
with the positions for x = 0.45). Since we do not know exactly all atomic positions at the 
compositions of interest, we cannot draw definite conclusions regarding phase separation. 



However, as pointed out earlier ||55|| , the relaxation of the lattice is very important and favors 
phase separation. 

The last comparison between the two types of superstructures we make is for x = 7/8. 
Since we do not include covalent corrections for Cu(l)-0(4) chains of intermediate length 
present in the HS (see Fig. 1), we drop E p i+E ch in this comparison. According to Table I, the 
magnitude of the neglected terms is ~ 3 eV. As in previous cases, only two-fold coordinated 
Cu ions are +1 and the rest are Cu +2 . As a consequence of the neglect of covalency, all O 
atoms of Cu02 planes are O -2 , and all chain 0(4) atoms are O - , except one of the four 
0(4) atoms of the 2\/2 x 2\/2 unit cell, nearest to the additional 0(4) vacancy, which is 0~ 2 
(we have chosen it inside the short chain). For the CS, at two-fold coordinated Cu(l) + 
is -9.14 V and that at four-fold coordinated Cu(l) +2 lies near -19.5 V. The potential at 
all apical 0(1) (chain 0(4)) ions lie near 22 V (16.7 V), with little variation with distance 
to the Cu(l)-0 vacancy chains. For HS, $ at four-fold coordinated Cu(l) +2 ions is -23.10 
V, except at the one nearest-neighbor to the 0(4)~ 2 , which amounts to -27.48 V. At the 
three-fold coordinated Cu(l) +2 ion nearest-neighbor to the 0(4) -2 ion, $ is also -23.10 V, 
and at the other three-fold coordinated Cu(l) +2 ion of the unit cell, = —18.73 V. At the 
0(4) -2 ions is $ = 17.57 V. The potential at the other 0(4) sites (occupied by O -1 ions) 
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vary between 12.42 V and 14.57 V. 

The energy of the HS for x = 7/8 is E = Eu a d — —292.44 eV. This is less than the 
energy E = Euad + -Ea for the CS, even at the largest possible A = Amux = 38.29 eV, 
for which E = —291.40 eV. Including E p i + E c h this energy decreases to -292.92 eV, but in 
principle one expects a similar decrease for the HS. Also for A = A min = 31.34 eV, even 
including E p i + E c h, the energy of the CS is -291.85 eV, larger than that of the HS. We 
conclude that our model supports the latter as the ground-state superstructure. 



IV. DISCUSSION. 

We have studied the interplay between the electronic and atomic structures of 
YBa 2 Cu30 6+x by a novel approach based on the Madelung energy and the cost of the 
charge transfer process Cu + + O — > Cu +2 + O -2 . The effect of covalency is included 
as a correction. This approach is motivated by the fact that while first-principles calcula- 
tions fail to describe the semiconducting systems, and the correlation energy they neglect 
is near 0.6 eV and depends on the oxygen ordering [|TJ, the strong-coupling models used so 
far ]l], 22| , 2J - gg] neglect long-range repulsions and depend on parameters which are not well 



known. A particular difficulty of these Hubbard-type models when applied to defects or 
systems with low symmetry (as the 1 x 8 or 2\f2 x 2\/2 unit cells) is that they require a 
large number of uncertain parameters (on-site energies at sites non-equivalent by symme- 
try for example) to describe the problem accurately enough. It is also difficult to include 
long-range repulsions in exact Lanczos diagonalizations of Hubbard like models. This was 
done by Riera and Dagotto for a generalized three-band Hubbard model in one and two 
dimensions ||56|| . However the two-dimensional results are incorrect due to subtleties in the 
use of boundary conditions, and in general different cluster sizes and shapes should be used 

We addressed the issue of the stability of "chain structures" (CS) of unit cell 1 x n 
(n integer) in comparison with 2 a/2 x 2y/2 type of superstructures (HS). As mentioned in 
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Section I, for x < 0.4 and room temperatures, the experimental evidence is against CS. 
We mention here also photoconductivity experiments |58|-|6T|: illuminating semiconducting 
(non CS) films, the resistivity decreases strongly as a consequence of pumping holes to 
the superconducting Cu02 planes, and ordering in (presumably short) chains takes place 
jBtJ since these structures are energetically favored under the constraint of a sizeable hole 
occupancy in the planes (this can be inferred from the information on the different potentials 
given in the previous section, or the arguments given in Refs. HI1J6^|). When illumination 
ceases, the resistivity returns to the original high values in times characteristic of oxygen 
diffusion (see Fig. 2 of Ref. 15§H), showing that the true equilibrium state is not a CS. 
Nevertheless, it is still possible that at lower temperatures a phase transition takes place 
(difficult to detect because of the sluggish oxygen kinetics at low temperatures), and the 
ground state is a CS. The present results support this statement. It is reasonable to expect 
that HS are favored by entropy at moderate temperatures: for CS, the cost in energy for a 
displacement of an O atom to their nearest available positions, breaking the chains, is high, 
while this is not the case for HS pi] . In fact this entropy term is essential to explain the 
neutron-diffraction results for x ~ 3/8 ]2D| in terms of a 2 a/2 x 4v^2 superstructure |3T| . 



The structural model of Aligia, Garces and Bonadeo P| J3^J5T[ ] is based on Coulomb 
repulsions between any two basal-plane 0(4) ions, screened by free carriers and dielectric 
polarization. Except for high stabilization energy of the chains AE, CS are unstable within 
this model for x ~ 1/8 because of the large cost in 0-0 Coulomb energy required to arrange 
the 0(4) atoms in Cu(l)-0(4) chains. The present results show that the neglect of electronic 
screening in the semiconducting phase is incorrect: for x = 1/8 and CS, the resulting 0(4) 
charge (~ —1), is screened by their nearest-neighbor Cu(l) +2 ions, and the cost in Coulomb 
energy to build the CS is not so high. Instead, we obtain that the cost of putting a line of O 
vacancies in the x = 1 structure is not so efficiently screened and thus, for x ~ 1, CS are not 
favorable. Other results which should be revised concern the hole count in the planes and 
the 60K plateau in the superconducting critical temperature [fl,|IT|, ^ , [2"4) --|2"E|] . The present 
results suggest that the role of apical 0(1) atoms is not so important as previously assumed 
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|I| JTT|j24| , and that the positions of the atoms and relaxation of the lattice are crucial in the 
charge balance, and also in a possible separation in phases with different oxygen contents 



0- 

Dielectric screening and that of free carriers, neglected in the present approach, are likely 
to play an essential role in the problem. One of the attempts to include Madelung poten- 



tials in electronic calculations is that of Ohta, Tohyama and Maekawa [|5^| . The authors 
screened those potentials by the optical dielectric constant (assumed 3.5 for all systems) to 
obtain different parameters of a multiband Hubbard model for the superconducting Cu-based 
perovskites. Using these parameters, they have obtained charge-transfer gaps, exchange con- 
stants and other information in good agreement with experiment. One would be tempted 
to extend trivially this formalism to calculate the total energy of YBa 2 Cu306+a;, dividing all 
Madelung contributions by 3.5. This is clearly incorrect, since dielectric and metallic screen- 
ing should lower the total energy of the system (increasing its absolute value) with respect 
to the unscreened case, in spite of the fact that the magnitude of the effective interaction be- 
tween two defect charges at a distance large in comparison with the interatomic distance, is 
reduced by dielectric or metallic screening JTJ,|32|. A simple electrostatic calculation involving 
two charges surrounded by a small void sphere (to avoid divergences) in a dielectric medium, 
shows that the interaction of the defect charges with the immediate neighborhood causes the 
largest reduction of the total energy. Local distortions around added or vacant 0(4) atoms 
in YBa 2 Cu 3 6+x were calculated by Baetzold |k| and are significant. In the semiconduct- 
ing phase, for which the effect of free carriers can be neglected, a formalism which takes 



into account Madelung energies, atomic potentials and atomic polarizations exists |64| , and 
might be applied to YBa2Cus06+ x for x < 0.4, as an extension and improvement of the 
present approach. 
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TABLES 





Zi 


Pi (V) 


Y 


3 


-29.39 


Ba 


2 


-18.41 


Cu(2) 


2 


-26.22 


0(2) 


-1.80 


20.04 


0(3) 


-2 


21.73 


0(2) 


-2 


21.76 


0(3) 


-1.80 


20.01 


Cu(l) 


2 


-24.35 


0(1) 


-2 


20.11 


0(4) 


-1.39 


15.87 


E 


-291.32 


eV 


EMad 


-288.59 


eV 


E A 





eV 


E p i 


-1.24 


eV 


E c h 


-1.49 


eV 



Table I: Charges (Zj), potentials ($) at the different atomic sites, and different contri- 
butions to the total energy of YBa 2 Cu306+a; for x — 1, e = 2 eV and A < 46.05 eV. 
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Zi 


Pi (V) 


Cu(2) 


2 


-26.73 


0(2) 


-1 91 


20 34 


0(3) 


-2 


21.15 


0(2) 


-2 


21.16 


0(3) 


-1.91 


20.32 


Cu 2 (l) 


1 


-12.23 




2 


-24 Q4 


0q(1) 


-2 


20.18 


4 (1) 


-2 


20.14 


0(4) 


-1.37 


15.59 


E 


-293.86 


eV 


EMad 


-276.93 


eV 


E A 


-15.5 


eV 


E p i 


-0.71 


eV 


E c h 


-0.73 


eV 



Table II: Same as Table I for the CS of x = 1/2. The subscript of Cu(l) ions refer to its 
coordination, and that of apical 0(1) ions is the coordination of its nearest neighbor Cu(l). 
Parameters are e = 2 eV and A = 31 eV. 
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z 


di (V) 


Cu(2) 


2 


-28.05 


0(2) 


-2 


19.79 


Of 31 


_2 


19.81 


Cu(l) 


2 


-21.40 


0(1) 


-2 


21.78 


0(4) 


-2 


21.25 


E 


-294.03 


eV 


EMad 


-293.93 


eV 


E A 





eV 


E p i 





eV 


E c h 


-0.1 


eV 



Table III: Same as Table I for the HS ("herringbone" |27j) for x = 1/2, and A < 
40.29 - e. 
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Fig. 1. Oxygen ordered superstructures of the basal plane considered in this work, for 
different oxygen contents x. Left: CS, of unit cell lxn. Right: HS, of unit cell 2^2 x 2y/2. 
Crosses denote Cu(l) atoms and solid circles represent 0(4) atoms in the notation of Refs. 
H,[44|. 
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